QCD corrections to J/ip plus i] c production in e + e annihilation at y/s = 10.6 GeV 
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Next-to-Leading-Order(NLO) QCD corrections to J/ip plus T} c production in e + e~ annihilation at 
\fs = 10.6 GeV is calculated in this paper, and an analytic result is obtained. By choosing proper 
physical parameters, a K factor (ratio of NLO to LO) of about 2, which is in agreement with the 
result in Ref. is obtained. The plot of the K-factor vs the center-of-mass energy y/s shows that 
it is more difficult to obtain a convergent result from the perturbative QCD without resummation 
of ln(s/mc) terms as the %/s becomes larger. 
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I. INTRODUCTION 

Perturbative quantum chromodynamics calculations 
are essential in the effort to describe large momentum 
transfer processes. To apply it to heavy quarkonium 
physics, the nonrelativistic QCD (NRQCD) factoriza- 
tion approach[2j has been introduced. It allows con- 
sistent theoretical prediction to be made and to be im- 
proved perturbatively in the QCD coupling constant a s 
and the heavy-quark relative velocity v. However, the 
J/ip polarization measurement at Fermilab Tevatron in 
proton-antiproton collisions H and J/ip production in 
B-factories [J, H Q have shown that the leading or- 
der (LO) theoretical predictions in NRQCD could not 
match the experimental results. The large discrepancy 
was found in the double charm production in e + e~ anni- 
hilation at B factories. The exclusive production cross 
section of double charmonium in e + e~ — > J/ipr] c at 
yfs = 10.6 GeV measured by Belle @, i] is a[J/ip + 
Vc] x B n " [> 2} = (25.6 ± 2.8 ± 3.4) fb and by BABAR 
6} is a[J/ip + r] c ] x B^{> 2] = (17.6 ± 2.8^1) fb, where 
B Vc [> 2] denotes the branching fraction for the T] c de- 
caying into at least two charged tracks. Meanwhile, the 
NRQCD LO theoretical predictions in the QCD cou- 
pling constant a s and the charm-quark relative velocity 
v, given by Braaten and Lee Ml, Liu, He and Chao 3, 
and Hagiwara, Kou and Qiao [0] are about 2.3 ~ 5.5 fb, 
which is an order of magnitude smaller than the exper- 
imental results. Such a large discrepancy between ex- 
perimental results and theoretical predictions brings a 
challenge to the current understanding of charmonium 
production based on NRQCD. Many studies have been 
performed in order to resolve the problem. Braaten and 
Lee have shown that the relativistic corrections would 
increase the cross section by a factor of about 2, which 
boost the cross section to 7.4 fb. And the NLO QCD cor- 
rection of the process has been studied by Zhang, Gao 
and Chao [l|, which can enhance the cross section with 
a K factor (the ratio of NLO to LO) of about 2 and 
reduce the large discrepancy. Again the relativistic cor- 
rections have been studied by Bodwin, Kang, Kim, Lee 
and Yu [l(| and by He, Fan and Chao [ll|, which are 
significant, and when combined with the NLO QCD cor- 
rections, may resolve the large discrepancy. In Ref. [ll|, 



T(4s) — > J/ip + rj c was considered by Jia, but it's con- 
tribution is small. Ma and Si [13j treated the process by 
using light-cone method. A similar treatment was per- 
formed by Bondar and Chernyad [l4[ and Bodwin, Kang 
and Lee [15|. More detailed treatment, such as including 
the resummation of a class of relativistic correction, has 
been take into consideration by Bodwin and Lee and Yu 

M- 

Since the calculation of NLO QCD correction for this 
process is quite complicate and plays a very important 
role to explain the experimental data, in this paper we 
performed an independent calculation by u sing the pack- 
age Feynman Diagram Calculation (FDC) [lj with one- 
loop part built in and obtained analytic result. The nu- 
merical result is in agreement with the previous result in 
Ref. pj. 

This paper is organized as follows. In Sec. II, we give 
the LO cross section for the process. The calculation of 
NLO QCD corrections are described in Sec. III. In Sec. 
IV, numerical results are presented. The conclusion and 
discussion are given in Sec. V. In the Appendixes, some 
useful details are presented. 



II. THE LO CROSS SECTION 

There are four Feynman diagrams for this order, and 
two of which are shown in Fig.[TJ while the other two can 
be obtained by reversing the arrows of the quark lines. 
Momenta for the involved particles are labeled as 



e (pi) + e+(p 2 ) -> J/ip[Pa) + Vc{p4)- 



(1) 



In the nonrelativistic limit, we can use the NRQCD fac- 
torization formalism and obtain the square of the scat- 
tering amplitude as 
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FIG. 1: Feynman diagrams for LO. 

J>3 P,i Ps 




FIG. 2: All Feynman diagrams for NLO are divided into thir- 
teen groups, (a) includes the photon-quark vertex counter- 
term and corresponding loop diagrams; (b) and (c) are the 
gluon-quark vertex counter-term and corresponding loop dia- 
grams; (d) and (e) denote the counter-term and corresponding 
loop diagrams for the quark and gluon self-energy; (f) and (j) 
are diagrams that contain Coulomb singularity. Other dia- 
grams can be obtained by reversing the arrows of quark lines 
and/or change the places of J/ip and r) c . But notice that we 
can't change the places of J/tp and tj c in group (h) and (i). 



where e 

,J/il> 



^ is the electric charge of the charm quark. 

R's'^iO) and R^°{0) are the radial wave function at the 
origin of J/tp and r/ c . Notice that s in Eq. §5§ is used 
from now on. After the integration of phase space, the 
total cross section is 
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the renormalization of the QCD gauge coupling constant, 
the charm quark mass and field, and the gluon field. A 
similar renormalization scheme is chosen as in ref. [HI 
except that both light quarks and charm quark are in- 
cluded in the quark loop to obtain the renormalization 
constants. The renormalization constants of the charm 
quark mass Z m and field Z2, and the gluon field Z% are 
defined in the on-mass-shell(OS) scheme while that of 
the QCD gauge coupling Z g is defined in the modified- 
minimal-subtraction(MS) scheme: 
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where is Euler's constant, /3 = ~ ^F n f is 

the one-loop coefficient of the QCD beta function and 
n/ is the number of active quark flavors. There are 
three massless light quarks u,d,s, and one heavy quark 
c, so n/=4. In SU(3) C , color factors are given by 



, Cf 



§,CU = 3. And 0' o = O + (4/3)2> = 
(11/3)Ca — (i/3)Tpnif where nif = rif — 1 = 3 is the 
number of light quarks flavors. Actually in the NLO to- 
tal amplitude level, the terms proportion to 5Zz° s cancel 
each other, thus the result is independent of renormaliza- 
tion scheme of the gluon field. 

After having fixed our renormalization scheme and 
omitting diagrams that do not contribute, including 
counter-term diagrams, there are 80 NLO diagrams re- 
mained, which are shown in Fig. [2l It is divided into 13 
groups. Diagrams of group (/) and (j) that have a virtual 
gluon line connected with the quark pair in a meson lead 
to Coulomb singularity ~ 7r 2 /v, which can be isolated by 
introducing a small relative velocity v = \p c — p s \- The 
corresponding contribution is also of 0(a s ) and can be 
taken into the cc wave function renormalization! lfj as: 



III. THE NLO CROSS SECTION 

Since there is no 0(a s ) real process in NLO, we only 
need to calculate virtual corrections. Dimensional reg- 
ularization has been adopted for isolating the ultravio- 
let (UV) and infrared (IR) singularities. UV divergences 
from self-energy and triangle diagrams are canceled upon 
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A factor of 2 should be used since there are two bound 
states. After adding contributions from all the diagrams 
together, all the IR-divergent terms are canceled and the 
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total scattering amplitude is obtained as 
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with i^i(s) and if 2 (s) given by Eqs. (|A~l~j) and (|A~2l) . 

Meanwhile a s should be obtained from two-loop for- 
mula as 
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where ft = 346^/3 - 4(Cf + 5C A /3)T F n f is two-loop 
coefficient of the QCD beta function. From Eq. ([7]) the 
total cross section at NLO is: 
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IV. NUMERICAL RESULT 

Up to NLO, the value of the wave function at the origin 
of J/tp is related to the leptonic decay widths as: 



16 a s 
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|i? s J ^(0)| 2 , (10) 



and according to ref. @] , we can set R^ c (0) = R J S ^ (0) = 
R s (0). If we choose |i? s (0)| 2 = 0.978 GeV 3 and = 
0.338 GeV, then we get the numerical result as shown in 
TABLE. U which is consistent with the result in ref. [l|. 
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TABLE I: Cross sections with different charm quark mass m c 
and renormalization scale /i. y^so = 10.6 GeV is the center- 
of-mass energy. 
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FIG. 3: Cross sections as function of the renormalization scale 
fi with |-R S (0)| 2 = 0.978 GeV 3 , A = 0.338 GeV and center-of- 
mass energy 10.6 GeV. The charm quark mass is chosen as 
1.4 GeV (upper curves) and 1.5 GeV (lower curves). 
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FIG. 4: Cross sections as function of the center-of-mass en- 
ergy with |i? s (0)| 2 = 0.978 GeV 3 and A = 0.338 GeV. Renor- 
malization scale fi is set at half of the center-of-mass energy 
and m c = 1.5 GeV. 



V. CONCLUSION 

We calculated the NLO QCD correction of J/ i/j plus rj c 
production in e + e~ annihilation at center-of-mass energy 
10.6 GeV. The method of dimensional regularization is 
taken to deal with the UV and IR singularities, and the 



Coulomb singularity is isolated by a small relative veloc- 
ity v between the charm quark pair in the meson and 
absorbed into the cc bound state wave function. After 
taking all one-loop diagrams into account, an analytic 
finite result is obtained. By choosing proper physical 
parameters, we get a K factor (ratio of NLO to LO) of 
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FIG. 5: K factor as function of the center-of-mass energy with 
|.R S (0)| 2 = 0.978 GeV 3 and A = 0.338 GeV. Renormalization 
scale is set at half of the center-of-mass energy and m c = 
1.5 GeV. 



about 2, which is consistent with ref. p|. It decreases 
the great discrepancy between theory and experiment. 
From Fig. [3l it could be found that the dependence on 
the renormalization scale \x has not been improved in 
NLO calculation. The plot of the total cross section vs 
the center-of-mass energy of e + e~ in Fig. |4] behaves as 
expected. But the plot of the K-factor vs the center- 
of-mass energy of e + e~ in Fig. [5] shows that it is more 
difficult to obtain the convergent result from the pertur- 

bative QCD without resummation of In — » terms as the 

mj. 

center-of-mass energy of e + e~ becomes larger. 
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APPENDIX A: THE DEFINITION OF Ki AND K 2 

In this section, the definition K\(s) and K2 (s) used in 
the Eq. ([7]) are presented. 



-24 - 79s - 68s 2 s 2 - 3s + 16 3s -34s 3 + 193s 2 - 342s + 160 

+ h — 7-77 7\ J 2 777 7v7 + J3 — 



3s(2s + 1) 



2(s-4) J 2(s-2) 



4(s-4)(s-2) 



-41s 2 - 194s + 64 8s 2 - 21s - 8 65s 2 - 302s + 64 -3s 2 - 4s 
+h 32(s-4) + h 2(s-4) + h 32(s-4) + h 2(s - 4) + &1 



4zi(8- 7s) 
s 2 (s-4) 



-64s 4 + 406s 3 + lis 2 - 335s - 120 



(s-4)(2s + l) 2 
-32s 4 - 110s 3 + 639s 2 + 696s + 172 



«6 



2^(8 -7s) z 2 (-130 + 19s + 8s 2 ) 



s 2 (s-4) 



s(s - 4) 



(s-4)(2s + l) 2 



a 7 - 



4z x (7s - 8) 2z 2 (-8s 2 - 19s + 130) 



s 2 (s-4) 



a 9 - 



8(8-4) 



(Al) 



-7s + 8 -8s 2 -19s + 130 32s 4 + 110s 3 - 639s 2 - 696s - 172 (65s 2 - 302s + 64) 
K 2 (s) = 2zi— 77-7 +z 2 — 1 777777 ; lN9 h Z\ 7-7- — a i3 



(s - 4)s 2 



(s - 4)s 



(s-4)(2s + l) 2 



4s(s-4) 2 



(21s 3 - 64s 2 - 1396s + 2688) (73s 3 - 96s 2 + 444s - 896) 
777 777-7 77T5 ai + Z\ — —— — 03 



(s - 2)s(s - 4) 2 



2)s(s-4) 2 



(193s 3 - 1292s 2 + 2548s - 1280) (79s 2 - 82s + 64) (8s 2 - 21s - 8) 

" 2 1 777 777 — 7 7TT5 a 5 + Z\ —. 7TT5 06 - 4z\ -, 7775 a 7 



-Z\- 



8(s-2)s(s-4) 2 
(34s 3 - 193s 2 + 342s - 160) 



(aio - a s ) + z\ 



8s(s-4) 2 u 1 s(s-4) 2 
(183s 3 - 1624s 2 + 3316s - 1408) 



2(s-2)s(s-4) 2 v 1U OJ 1 8(s-2)s(s-4) 2 
Where all the variables used in Ki(s) and i^ 2 (s) are defined as: 
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APPENDIX B: THE RESULTS FOR ALL THE 
SCALAR INTEGRALS 

In this section, we present the results of scalar inte- 
grals. Functions A to E denote one to five point scalar 
integrals and variables of the functions are written as 
T(pQ,mo, . . . ,p n ,m n ) where pi and m, denote parame- 
ters of the nth propagator Ni — (q + Pi) 2 — mf + it. q is 



the loop momentum and a factor of 
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is taken away from all the scalar integrals. 

We have developed a full series of methods in calculat- 
ing tensor and scalar integrals with dimensional regular- 
ization and realized it in FDC[l7|. A paper about these 
methods are in preparation. All the scalar integrals are 
calculated analytically by using FDC and the results are 
shown in the following. 
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And here is the results for Fi where /, and a, are defined 
as before. 
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[2a x (s + 1) + 2(-a 6 s + i?rs + 2s + 1)] 



i(2a 1 + l) 



1 



F] 
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2 [8tti(-s + z 2 ) + 8z 2 (2a 9 - a 6 ) + 8s(a 6 + 1)] 
(-32ai + 16) 



F 17 = -=-[167ri(— s + z 2 ) + 16z 2 (2ag — a 6 ) 
s^ 

+16(a6S + 2ais)] 
1 



F 



is 



32?ri(s - 2)(s - z 2 ) + 32z 2 (s - 2)(-2a 9 + a 6 ) 



-32a 6 s(s - 2) - 64ais 2 + 64s] 



(B7) 



APPENDIX C: THE SCHEME TO TREAT 7s IN 
D-DIMENSION 



As we all know, there is no explicit definition for 75 
in D dimensions. Usually the following relations is used 
when one encounters 75 in D dimensions: 

{75,7mI = 0, (CI) 

Tr(78plp2P3P4) = ^^af3PiP2P3P4- ( C2 ) 

Notice that en Va p goes zero when any of it's indices is 
out of the 4 dimensions. 



While calculating the trace of the product of sev- 
eral matrices that contain 75 and indices in D dimen- 
sions, different ways may l ead t o different results. For 




FIG. 6: Trace calculation including 75. 



example, when calculating the trace of matrix M = 
757uPiP2PkP47^) we have two different routes as shown 
in Fig. [SJ 

1. In route A±, the summing up of the index fi does 
not go across 75. 



Tr(M)| Al = Tr 75 



-2p 3 pip 2 P4 + 2p4Pip 2 P3 



2. In route A 2 , it does go across 75. 



(C3) 



Tr(M)| A2 = -Tr 757 M 7 M pip 2 p3p 4 



DTr yy 5 PiP2P3P4 



-MDe^app^p^PsP^ 



(C4) 



It is easy to find that Tr(M)\ Al = Tr(M)\ A2 in 4 dimen- 
sions, but in D dimensions they're different to each other. 
So we should always take same route when dealing with 
traces containing 75 in order to keep our final finite result 
consistent. In FDC, route A\ is taken. 
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